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ABSTRACT 



The most massive Galactic globular cluster, lo Centauri, appears to have mul- 
tiple populations. Its bluest main sequence and extreme horizontal branch stars are 
suggested to have the common origin, that is, an extremely high helium abundance 
of y ~ 0.4. The high helium abundance is most often attributed to asymptotic gi- 
ant branch (AGB) stars. In this study we test the AGB hypothesis. We simulate the 
maximum-AGB models where the impact of AGB stars is maximised by assuming 
that supernova explosions do not affect the chemical evolution of the proto cloud. 
We compare the enrichment history of helium, metals, carbon and nitrogen to the 
observed values. Even under the most generous condition, the maximum-AGB models 
fail to reproduce the large values of helium Y ^ 0.4 and helium enrichment parame- 
ter AY/AZ ^ 70 which were deduced from the colour-magnitude diagram fits. They 
also fail to reproduce the C and N contents of the blue population spectroscopically 
determined. We conclude that the AGB scenario with the canonical stellar evolution 
theory cannot explain the observational constraints and that the self chemical enrich- 
ment does not provide a viable solution. Alternative processes are desperately called 
for. 

Key words: galaxies: globular clusters — individual (w Centauri) 



1 INTRODUCTION 



The origin of the blue population in the most massive 
Galactic globular cluster (GC), lj Centauri, is at the cen- 
tre of debat e after the dis c overy of i ts multiple popu- 
lation s (e.g. lAndersonI 1 19971 : iLee et all 1 19991 : iBedin et all 
[2003). The explanation of the blue main sequence (bMS) 
of OJ Cen by a huge excess of the primordial helium abun- 
dance. Ay ~ 1.2, was implied via the colour-magnitude di- 
agram (CMD) fitting (Norris 2004). This bMS, which con- 
tains ~ 30 percent of the cluster stars, is roughly 0.3 dex 
more metal-rich than the red main sequence (rMS) popu- 
lation (|Piotto et aLluOOSl ). Its blue colour despite its high 
metallicity is attributed to an extremely high value of he- 
lium. This explanation has further been reinforced by the 
fact that the same helium excess ca n reproduce the extended 
horizonta l bran ch (EHB) of w Cen l|D'Antona fc Caloi|[20o3 : 
iLee et al.l |2005|) . Follow-up spectroscopic and photometric 
inves tigations hav e been perform ed by many groups ( e.g . 
Kavser et al. 2005: ISoUima et al.ii2006: Stan ford et al. 20_07 



Villanova et all l2007l : Ivan Loon et al.l [20071 '). The extreme 



helium scenario appears to provide a successful solution 
to another clust er, NGC 2808, wh ich also shows multiple 
main sequences (jPiotto et al■ll2007^ and a pronounced EHB 



IjLee et ahl |2005| ). It has recently been claimed that the 
GCs with an EHB are among the most massive in the 
Milky Way, and the possible link between the mass of the 
cluster and the helium anomaly is heavily investigated on 
l|Recio-Blanco et al.l [JOO^ : iLee et al.l [2OO7I '). The existence 
of helium-rich GCs in the externa l galaxy, M87, has also 
been suggested (JKavirai et al.ll20(yn ). implying that this phe- 
nomenon may be universal. 

The extreme helium scenario, however, has been crit- 
icised. This is mainly because the helium-enrichment pa- 
rameter AY/AZ corresponding to proposed helium abun- 
dance r-^ 70 or even higher, which is more than an or- 
der of magnitude larger than obse rvational and theoreti- 
cal value s, AY/AZ = 1 -5 (e.g. |Fe rnandcs et all Il996l : 
Hagcl & Portinaril Il998l : I Jimenez etal . 2003). For possi- 
ble origins for the extreme helium abundance, asymp- 
totic giant branch (AGB) stars, massive stars, and 
Type I I supe r novae (SN II) h a ve be e n widely discussed 
llNorrid |2004'; DAntona et al. '2005'; 'Maeder fc Mevnetl 
l2006f ). iBekki & Norris (2006) however demonstrated that 
such candidates cannot produce the amount of helium re- 
quired, for ordinary initial mass functions (IMFs) within 
the scheme of a closed-box self enrichment. More recently, 
IChoi fc Yil (120071 ) showed that essentially no population can 
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produce such a high value of Ay/AZ via self-enrichment 
processes regardless of the shape of IMF. 

Supernovae are not helpful for generating a high value 
of Ay/AZ because they produce metals as well as helium. 
The ejecta from AGB stars is generally believe d to show the 
large s t helium enrichm ent parameters (e.g., see llzzard et al.l 
I2OO4I : IChoi fc YilboOTi ). However, the AGB hypothesis is 
found to be unsuccessful for matching both the high helium 
abundance and the CNO properties constrained from obser- 
vations eve n when the supernova effect is removed from the 
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calculation IjKarakas et al.ll2006l : iBekki et al.l 
W e further ela b orate on the works of 
120061 '): iBekki et al] (|2007l ') by testing the maximum- AGB 
hypothesis where the impact of AGB stars is maximised by 
assuming that massive star ejecta of unknown mass range 
escape the small potential wells of proto clouds. In this let- 
ter we search for the parameter space that satisfies all the 
observational constraints including spectroscopic measure- 
ments. 



2 THE CHEMICAL EVOLUTION MODELS 

2.1 Introduction 

We adopt the chemical enrichment cod e described i n 
IChoi fc Yil (|2007l ) following the formahsm of lTinslevI (|l980l ). 
The two-component model assuming instant mixing and 
cooling traces the net metallicity Z, helium Y , carbon and 
nitrogen contents of the gas and stellar components. 

The stellar mass, Ms{t), and the cold gas mass, Mg{t), 
are normalised to the initial system mass, 

Msit) 



^J.s{t) 



MaW 



MM 



(1) 



MtotiO)' '^"'^' ~ MtotiO)' 

In this study, we assume that the formation of the for- 
mer (rMS) population is essentially instant: 

^{t / 0) = 0. (2) 

The evolution of the gas mass is given by 

^ = E{t)~m (3) 

where the ejecta gas mass at time t, E{t), is defined as 
E(t) = I dmcj>{m){m~Wm)tp{t-Tm). (4) 



where (pi^ denotes the initial mass function (IMF), Wm 
the remnant mass for a star with main sequence mass m. 



Tm the lifetime of a star of mass m, and niu 



the upper 



mass cut in the initial mass function. The Scalo IMF IjScald 
1 198a ) with cutoffs at 0.1 and IOOMq is assumed, and the 
remn ant mass Wm and the l ifetim e of a star r™ are adopted 
from iFerreras fc SiJ^ (I2OO0I . I2OOII ) . 

The equation for the evolution of metallicity in gas is 
given by 

d(ZgfJ.g) 



dt 



-tl){t) + Ez{t) 



(5) 



where the mass of ejected metal at time f, Ez{t), is defined 
as 

Ez{t) = I dm(j){m)[mpmtp{t - Tm) + {m — Wm)ip{t ~ 

Jmt 



where pm denotes the newly synthesized and ejected metal 
(or helium, C, N) mass fraction to initial stellar mass. We 
will discuss this in greater detail in §2.2. 

The initial chemical properties of the models are 
adopted from the observational values for the rMS popu- 
lation of uj Cen. The initial metallicity Zq = 0.001 and the 
initial C, O abundance [C/M]= .0, [N/M]= 1 . are ft-om 
the spectroscopic measurements of IPiotto et al.l (2005), and 
the initial h elium abu n dance Yq = 0.232 is deduced from the 
CMD fit in lLeeet all l|2005h . 

2.2 Chemical Yields 

The chemical yields (pm), defined as the mass fraction of 
a star of mass m that is newly converted to metals or 
helium and ejected, are the most important input param- 
eter to the chemical evolution of a population. We basi- 
ca lly adop t the helium, carbon and metal pm predictions 
of lMaedej (| 19921 ) for the mass range 9- WOMp,. For nitro - 
gen, we adopt the pm prediction of iNomoto et al.l ([1997]). 
For Mq ^ 6 we adopt th e heli um, carbon, nitrogen, and 
metal yields from iHerw ig] KOOJ). We confirm that the use 
of different yields (e .g. Ivan den Hoek fc Groenewegenlll997l : 
IVentura et al.ll2002l ) does not make any notable difference to 
our conclusion. We use the chemical yields from the hermi- 
tian fits to the pm predictions. In Figure 1, we show the the- 
oretical ejecta abundance predicted by our fitting functions 
from the metal-poor {Z = 0.001) stars of mass 0.1 — IOOA/0 . 
The chemical abundances predicted and the resulting he- 
lium enrichment parameters are presented. As shown in the 
middle panel, the highest peak around 6 Mq in the helium 
enrichment parameter diagram reaches the value proposed 
for the bMS population. Indeed, the AGB ejecta appears to 
be the most promising candidate. Supernovae on the other 
hand produce more metals than helium and hence are not 
helpful for generating a high helium-to-metal population. 



2.3 The Maximum- AGB model 

It has already been shown that a closed-box system can- 
not achieve such high values of helium and helium enrich- 
ment parameter a s implied for the blue population of uj Cen 
IChoi fc Yill2007l '). This is mainly because of the ill effect of 
supernovae on AY/AZ. 

In realistic models for small potential wells, the rem- 
nants of sup ernovae are believed to escape the system 
(|Larsonlll974l ). Such explosions would probably remove the 
remnant gas in the potential as well after the first star for- 
mation episode. But in general, the escape of the supernova 
materials is likely to be easier than the escape of the rem- 
nant gas. So we set up ad hoc models in which the super- 
nova ejecta escape the system without affecting the chemical 
properties of the remnant gas, hence maximising the chemi- 
cal influence of AGB stars. We call this "the maximum- AGB 
models" . 

Our prescription has two input parameters. The first 
free parameter, niesc, indicates the critical mass above which 
the stellar mass ejecta escape the potential well without af- 
fecting the chemical abundance of the remaining gas. This 
parameter is likely a function of the mass of the system be- 
^a^se tjiCjpotpjD^tial well determines whether the supernova- 
escape the gravitational potential of the 
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Figure 1. The theoretical eje cta a bundance calc ulated from the 
stellar yields of lMaedeJ l ll992D and lHerwij ||2004|) . Top: The he- 
lium and metal abundances of the ejecta from the metal-poor 
{Z = 0.001) stars of mass 0.1 - WOMq. Middle: The resulting 
helium enrichment parameters. Bottom: The carbon and Nitrogen 
abundances of the ejecta from the stars wit h the initial compo - 
sitions [C/M] = 0.0, [N/M] = 1.0 following [Piotto et al.l ^0^. 



system or not. iKarakas et al.l l|2006l ) assumed that the sys- 
tem can retain the ejecta from the stars with m ^ 6.5Mq. 
We explore an expanded parameter space for 5Mq ^ rrissc ^ 
99Mq . The scenario of Karakas et al. is equivalent to our sce- 
nario with TUesc = 6.5Mq. The model with the largest ruesc 
represents the ordinary chemical evolution model which 
takes account of all the mass ejecta from stars when cal- 
culating the metallicity of the next generations. The model 
with the smallest rriesc shows the extreme condition that the 
ejecta from AGB stars dominates the chemical properties 
of the next generations, hence the maximum-AGB model. 
Then we define E{t), the ejecta gas mass at time i as 
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Figure 2. The chemical abundances of the second stellar genera- 
tion resulting from the maximum-AGB chemical evolution model. 
The results are shown for different values of initial starburst mass 
fractions, RrMS = 0.5-1.0 with AR^ms = 0.05. Dashed/solid 
lines show the models that are incompatible/compatible with the 
empirical mass fraction of the blue population (30 percent), for 
the age difference of 1 Gyr. The shaded region shows the results 
for the age difference of 3 Gyr. The chemical abundances of the 
blue population of uj Cen are shown by the black dotted hori- 
zontal line in each panel. Top: The helium enrichment parame- 
ter. The mass fractions of the blue population are shown for the 
models that are incompatible with the empirical requirement (30 
percent). Middle: The helium content. Bottom: The metallicity. 



of stars is solely influenced by the ejecta from the former 
generation. Then we define the star formation rate as 



^(i) = RrMS S{t) 



(9) 



E(t) = I dm4>{m){m - Wm)tp{t - Tm) , (7) 

J mt 

and Ez{t), the mass of the ejected metal at time t as 

Ez{t) = XT"" dm(j){m) [mp,„V(i - Tm) + (m - Wm)t{t - r„) 

Another free parameter is RrMS which indicates the 
mass ratio of the initial starburst i.e., for the former, red 
population of uj Cen to the total system mass. We assume 
that the formation of the former (red) population is essen- 
tially instant and regulate the fraction of the initial star 
formation by using the free parameter RrMS- We calculate 
models for 0.5 < Rtms ^ 1-0. The model with RrMS = 1-0 
means all gas goes into the star formation at time t — and 
the metallicity of the gas that forms the next generation 



where 5{t) is a Dirac delta function. 

In effect we are testing a self-enrichment hypothesis 
where the helium enrichment of the blue population is solely 
due to the ejecta from the former (red) population. We 
calculate the gas enrichment history for the predicted age 
V differencJ; bjOtwe en the tw o extreme populations of lu Cen, 
''i.e. 1-3 Gyr l|Lee et al . 2005; Stanford ct al. 2006). We 
also consider the observed number fraction of the bMS stars 
(~ SOpercent) as an additional constraint. The rMS-to-bMS 
ratio, 7:3, may not be a strong constraint if the bMS and rMS 
populations formed in geographically different regions. For 
example, the bMS stars are obser ved to be more cen trally 
concentrated than the rMS stars (|Sollima et al.l 120071 ). and 
thus if one assumes that the original rMS stars had a bet- 
ter chance of escaping from the system potential during the 
dynamical encounters with the Milky Way galaxy over the 
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Figure 3. The same as Figure 2, but for the spectroscopic mea- 
surements [C/M] (top) and [N/M] (bottom). 



Hubble time, the original mass ratio may have been higher. 
We discuss the impact of this uncertainty in §3. 



3 RESULTS 

Figure 2 shows the helium and metal contents and the he- 
lium enrichment parameter of the remaining gas after a 1- 
Gyr evolution. The chemical properties of the models are 
presented for our m^ac - RtMS parameter space. The mod- 
els assuming Ai ~ 1 Gyr (that is, the age difference be- 
tween the old red and the young blue populations) are 
shown by lines, while the shaded region shows the results 
for At ~ 3 Gyr which is close to the upp er limit on the 
age difference between the two populations IjStanford et al.l 
12006 ) . The solid/dashed lines show the models that are com- 
patible/incompatible with the observed mass ratio, 7:3, re- 
spectively. As shown in the top panel of Fig 2, the model 
with rriesc = TM© and RtMS = 1.0 has the largest value 
of /SY//S.Z ~ 10.5. This val ue is much higher than the 
canonical value, AY/AZ ~ 2 (|Pagel fc Portinaril [l998l ') yet 
still much lower than the value we are seeking for (i.e., 
Ay/AZ ~ 70). In addition, the maximum helium con- 
tent which can be achieved in this model is Y ~ 0.33 (for 
Tiesc 5S 5OM0 and RtMS = 1-0), also much smaller than 
Y ~ 0.4 estimated by the CMD analysis of to Cen. Even 
when we assume an extreme condition, rriEsc = 5 Mq, that is 
hardly plausible in terms of physics, the extremely high val- 
ues of helium and helium enrichment parameter suggested 
for the bMS cannot be reproduced. 

As mentioned in §2, the mass fraction of the bMS popu- 
lation (30 percent) may be an upper limit. As seen in the top 
and middle panels of Figure 2, models cannot reach the re- 
quired values even when we lower the bMS mass fraction 
constraint by a significant factor. The situation becomes 
even worse when we compare the models to the spectro- 
scopic line strengths. In Figure 3, we present the resulting 
carbon and nitrogen contents of the models for the same 
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Figure 4. The calculated chemical abundance of the ejecta for 
the entire (RrMS = 0-5 ~ l-Oi ^^ = 0.05) parameter spaces of 
the maximum-AGB model is represented. Filled circle and cross 
represents the resulting value of ruesc = 5,15, 99Mq model. Filled 
circle indicates that the condition of mass ratio between rMS and 
bMS is satisfied while cross denotes the model which do not meet 
this condition. Shaded region shows the viable region of the total 
parameter space, rriesc = 5 — 99Mq, which satisfies the mass 
ratio condition. Reference values of blue population of uj Cen are 
marked(filled star) in each panel. 



niesc and Rtms parameter space as in Figure 2. This result 
is important since the reference chemical properties of the 
bMS are directly measured via spectroscopy, while the he- 
lium abundance and the helium enrichment parameter are 
deduced from CMD fits. The models with mesc ~ 7Mq show 
the largest value of AY/AZ (Figure 2 top panel) but do not 
match the observed carbon content of the bMS (Figure 3 
top panel). The horizonta l, shaded band show s the observed 
value for the bMS stars l|Piotto et al.l '20051. However, the 
errors are not given and ca n be more significant than we esti- 
mate d from the literature (jPiotto et al.ll2005l : [Stanford et al.l 
[2007]). I n case of the red p opulation, [N/M] is poorly con- 
strained; |Piotto_et_aL| (120051 ) mentioned that values of [N/M] 
^1.0 could also be compatible. We set [N/M] = 1 for the 
red population. This assumption on the nitrogen abundance 
of the rMS population may appear unphysically high con- 
sidering the typical nitrogen production of stars, that is, an 
order of [N/M]--^ —1 through 0, as shown in Figure 1. The 
use of lower values of [N/M]rMS results in a substantially 
greater mismatch in [N/M]bAfS in this figure. 

We present the resulting carbon, nitrogen and helium 
abundances simultaneously in Figure 4. The three cases of 
rriesc ~ 5, 15, 99Mq are shown. Filled circle (cross) show the 
models that are consistent (inconsistent) with the mass ratio 
constraint, respectively. The shaded region shows the the 
parameter space that satisfies the mass ratio constraint. The 
observed chemical properties of the blue population stars of 
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Figure 5. The same as Figure 4 but for tlie [C/M] and [N/M] 
space. 



OJ Cen (filled star in each panel) are not reproduced under 
any circumstance. Note the logarithmic scale of the axes. 

Figure 5 sho ws that the models fail to match the spec- 
troscopic data of IPiotto et al.l (|2005l ). The empirical value 
(star symbol with errors) measured using the VLT is in rea- 
sonable agreement w ith the more recent measurements of 
[Stanford et al.l (|2007l ). As mentioned, the models are based 
on the assumption of [C/M] = 0.0 and [N/M] = 1.0 for the 
red population. The model values of [N/M]f,MS may be con- 
sidered upper limits. The mismatch is significant. 



still much lower than suggested by the CMD fits. Even if 
we admit that the estimation for the helium abundance 
via CMD fits is extremely difficult and so uncertain, the 
maximum- AGB scenario cannot reproduce the carbon and 
nitrogen properties spectroscopically observed and hence 
probably more robust, either. 

For an al ternative solution, surface pollu- 
tion scenarios (ITsuiimoto. Shigevama. fc Sudal 120071 : 



iNewsham fc Terndrod [2007|) have been suggested and 
deserve attention especially because they demand only 
mildly-enhanced helium abund ance for t he bl ue popula- 
tion. From a different avenue, IChuzhovl (|2006l ) suggested 
that the primordial helium may have been concentrated 
in stellar-mass scales within minihalos of size roughly 
consistent with dwarf galaxies. This aspect is particularly 
appealing because uo Cen is often considered as a remnant 
of a former dwarf galaxy and the multiple main sequences 
are found in the most massive globular clusters in general. 
This scenario offers an attractive solution to the helium 
variation in first stars. However, it remains unclear how 
long the extreme-helium clumps could last. The blue 
population of uj Cen appears to be substantially enhanced 
in metals as well and hence hardly qualifies as first stars 
and y ounger than the re d population by a couple of billion 
years. I Choi fc Yil (|2007l ) proposed the fluctuation in the 
chemical properties of the initial starburst clumps during 
and after the first star formation epoch. Given the impact 
of the issue to the understanding of the galaxy formation, 
a more detailed and realistic model calculation is necessary 
until such scenarios become convincing. 



4 DISCUSSION 

We have used a chemical enrichment model to demonstrate 
that the AGB ejecta of the red population cannot have pro- 
duced the blue population of ui Cen. We tested a hypothesis 
in which the massive star ejecta escape the potential well so 
that the AGB effect can be maximised. But even in this ad 
hoc scenario, models miserably fail to match the chemical 
properties of the blue population of uj Cen. We use the bMS 
mass fraction (30 percent) and the age difference between 
the two populations (1-3 Gyr) as constraints, but changing 
them does not make the match any bette r. We as a result 
confirm the result of iKarakas et al.l l|2006l '). 

The helium content of the AGB ejecta is basically too 
small to be consistent with the value suggested for the 
blue population o f lo Cen, Y'-^ . 4, as we discussed earlier 
l|Choi fc Yill2007l '). [Romano et al.1 l|2007h has also confirmed 
this using various yields from indep e ndent studies (e.g. 
van den Hoek fc Groenewegen 1997 : Mevnet fc Maeden 



20021 : iHirschil l2006l : iMevnet. Ekstrom. fc MaedeJ 120061 ) 
Such a high level of h elium excess can be achieved only 
through massive stars (iNorrisI |2004| : iD'Antona et al.1 l2005l : 
iMaeder fc Mevnetll200q ). However, such massive stars also 
produce a large amoun t of metals and thus resulting in a 
small value of AY/A^ (|Choi fc Yill2007l '). In this respect, a 
variation in the initial mass function is useless. The escape 
of the massive star ejecta through supernova-driven winds 
would lead to a larger AY/AZ roughly up to 10; yet. 
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